Application of relativistic mean field and effective field theory densities to scattering observables for 

Ca isotopes 
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In the frame work of relativistic mean field (RMF) theory, we have calculated the density distribution of 
protons and neutrons for 4 °. 42 > 44 > 48 Q a with NL3 and G2 parameter sets. The microscopic proton-nucleus 
optical potentials forp _|_ 4 °. 42 > 44 . 48 Q a systems are evaluated from the Dirac NN-scattering amplitude and the 
density of the target nucleus using Relativistic-Love-Franey and McNeil-Ray-Wallace parametrizations. We 
have estimated the scattering observables, such as elastic differential scattering cross-section, analyzing power 
and the spin observables with the relativistic impulse approximation (RIA). The results have been compared 
with the experimental data for few selective cases and find that the use of density as well as the scattering matrix 
parametrizations are crucial for the theoretical prediction. 

PACS numbers: 25.60.Bx, 25.10.+S, 26.60.Gj, 26.30.Ef,26.30.Hj 



I. INTRODUCTION 

Study of the nuclear reactions is a challenging subject of 
nuclear physics both in theory and laboratory. This is useful to 
explain the nuclear structure of stable as well as exotic nuclei. 
The Nucleon- Nucleus interaction provides a wide source of 
information to determine the nuclear structure including spin, 
isospin, momenta, densities and gives a clear path towards the 
formation of exotic nuclei in the laboratory. In this context the 
study of elastic scattering of Nucleon-Nucleus is more inter- 
esting than that of Nucleus-Nucleus at different energies. One 
of the theoretical method to study such type of reactions is 
the "Relativistic Impulse Approximation" (RIA). It is a micro- 
scopic theory where the Dirac optical potential is constructed 
from the Lorentz invariant Nucleon-Nucleon (NN) amplitudes 
obtained from relativistic meson exchange models. The basic 
ingredients in this approach are the NN scattering amplitude 
and the nuclear scalar and vector densities [ 1 ] of the target nu- 
cleus. This approach can be extended to elastic scattering of 
composite particles [[2]]. In this context proton-Nucleus (p-A) 
elastic scattering is of particular interest because of its relative 
simplicity with which it provides a satisfactory description of 
the reaction dynamics. 

One useful application of RIA is to generate microscopic 
optical potential to study the elastic and inelastic scattering of 
nucleons for unstable proton- and/or neutron- rich nuclei. The 
RIA folding procedure can also be extended to calculate mi- 
croscopic optical potentials for exotic nuclei using relativistic 
mean field formalism JUS]. 

The first theoretical introduction to elastic scattering was 
given by Chew [5] almost six decades ago. For a wide range 
of energy interval, Impulse Approximation (IA) produces the 
main qualitative description on quasi-elastic scattering for 
A < 64 nuclei |g|. During the same time, Glauber [7] studied 
the reaction dynamics of the composite system at low energies 
but this model is unable to predict extension of quasi-elastic 
scattering. Further the generalized Glauber formula and the 
unitarized impulse approximation [8, 9] were circumvented. 
But the development of RIA opens a path to study the above 



mentioned scattering phenomenon for both the elastic and 
the quasi-elastic particles. This field of research was further 
strengthened by the experimental evidences of cross-section 
and analyzing power for the scattering systems p + 12 C, 
p + 9 B and p + 16 O at 200 MeV, which were measured over a 
wide range of momentum transfer > 6fm~ 1 at IUCF [ lfj.ll ill . 
Recent study of proton nucleus elastic scattering within modi- 
fied (Coulomb) Glauber model fL2l[l3ll and global Dirac opti- 
cal potential [ 14] have motivated us to study the elastic scatter- 
ing phenomenon. For convenience, we consider Ca isotopes 
as targets and p as a projectile, because Ca satisfies the rela- 
tivistic mean field nuclear structure model accurately without 
recoil correction to the Dirac scattering equation. 

In the present paper, our aim is to calculate the nucleon- 
nucleus elastic differential scattering cross-section (4^) and 
other related physical quantities, like optical potential (U op t), 
analyzing power (A y ) and spin rotation parameter (Q— value) 
using relativistic mean field (RMF) and recently proposed ef- 
fective field theory motivated relativistic mean field (E-RMF) 
densities. These are obtained from the successful NL3 [15] 
and the advanced G2 lloll parameter sets, which are given 
in Section II. In the Sections III and IV, the details of tar- 
get densities folded with the NN-amplitude for various en- 
ergetic proton projectile with Relativistic-Love-Franey (RLF) 
1 17, [3 and McNeil-Ray-Wallace (MRW) parametrizations 
1 1911 for 4 °. 42 . 44 . 48 C' a are given. In these sections we have 
outlined the expressions for the differential elastic scattering 
cross-section, analyzing power and spin observables. Section 
V describes the results obtained from our calculations. Finally 
a brief summary and conclusions are given in the Section VI 
for the present work. 



II. THE RMF AND E-RMF FORMALISMS 

A documentation of RMF and E-RMF formalisms are avail- 
able in Refs. l20ll and flrjl I2TI1 respectively for both finite 
and infinite nuclear matter. Here only the energy density 
functional and associated expressions for the densities are 
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where the index a runs over all occupied states <p a (r) of the 
positive energy spectrum, $ = g s <fro(r), W = <7t,Vo(r), i? = 
g p b (r) and A = eA (r). 

The terms with <ju,, A, (3 S and /3„ take care of the effects 
related with the electromagnetic structure of the pion and the 
nucleon (see Ref. rf23to . Specifically, the constant g 1 con- 
cerns the coupling of the photon to the pions and the nucleons 
through the exchange of neutral vector mesons. The experi- 
mental value is g^/^n = 2.0. The constant A is needed to 
reproduce the magnetic moments of the nucleons, defined by 



(2) 



with A p = 1.793 and A„ = —1.913, the anomalous mag- 
netic moments of the proton and the neutron, respectively. 
The terms with f3 s and f3 v contribute to the charge radii of 
the nucleon ll23ll . 

The energy density contains tensor couplings, scalar-vector 
and vector-vector meson interactions in addition to the stan- 
dard scalar self-interactions K3 and K4. Thus, the E-RMF for- 
malism can be interpreted as a covariant formulation of den- 
sity functional theory as it contains all the higher order terms 
in the Lagrangian, obtained by expanding it in powers of the 
meson fields. The terms in the Lagrangian are kept finite by 
adjusting the parameters. Further insight into the concepts of 
the E-RMF model can be obtained from Ref. M23I1 . It may 
be noted that the standard RMF Lagrangian is obtained from 
that of the E-RMF by ignoring the vector-vector and scalar- 
vector cross interactions, and hence does not need a separate 
discussion. 

In each of the two formalisms (E-RMF and RMF), the set of 
coupled equations are solved numerically by a self-consistent 
iteration method. The baryon, scalar, isovector, proton and 
tensor densities are 
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These densities are obtained from the RMF and E-RMF for- 
malisms with NL3 d 1 511 and G2 B16I1 parametrizations. We 
refer the readers to Refs. I201 121 1 for numerical details and 
ground state equations for finite nuclei. 



in. THE NUCLEON-NUCLEON SCATTERING 
AMPLITUDE 

The non-linear relativistic impulse approximation (RIA) in- 
volves mainly two steps 124142611 of calculation. Basically a 
particular set of Lorentz covariant function lfl9ll . which multi- 
ply with the so called Fermi invariant Dirac matrix represent 
the Nucleon-Nucleon NN-scattering amplitudes. This func- 
tions are then folded with the target densities of proton and 
neutron from the relativistic Langragian for NL3 and G2 pa- 
rameter sets to produce a first order complex optical poten- 
tial. The invariant NN-scattering operator T can be written 
in terms of five complex functions (the five terms involves in 
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the proton-proton pp and neutron-neutron nn scattering). In 
general RIA, the function T can be expressed as fl24 - l26ll . 



PS 



A (i) 



(9) 



where X L stands for Dirac operator and (0) and (1) for the in- 
cident and struck nucleons respectively. The S, V, T, A and 
PS stands for scalar, vector, axial vector, tensor and pseu- 
doscalar. The dot product (.) implies all Lorentz indices are 
contracted. The Dirac spinor is defined the initial and final two 
nucleons by taking the matrix elements of T, which represent 
the NN-scattering amplitudes. The function T L are deter- 
mined by equating the resultant amplitude (in center of mass 
frame) to the empirical amplitude, which is conventionally ex- 
pressed in term of the non-relativistic Wolfenstein amplitudes 
A\ , A2, ....A5 lfl9ll . Since there are five complex invariant am- 
plitudes and Ai, A%, ....A$ are five Wolfenstein amplitudes, 
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the relation among them is determined by a 5 x 5 non-singular 
matrix, whose inversion is straight forward. However T is an 
operator in the two particle Dirac space and the component 
cancelled out due to isospin and parity invariance and we get 
only 44 components [27]. From the above, it is clear, to spec- 
ify that T is not unique. In other words, there are infinite 
number of operators T with same five on-shell but different 
negative (energy) elements. The expression of T cannot pre- 
dict reasonanle result at lower energy region. To avoid the lim- 
itation, the pseudoscalar T PS is replaced by the pseudovector 
invariant, and is expressed as, 
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The meson-nucleon couplings are complex, with a real part 
gf and an imaginary part gf, which can be decomposed into 
two parts, 
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where t(E) is the lowest order meson and T is the total isospin 
of the two nucleon state. The calculation of the one-meson- 
exchange from Feynman diagram IU7I1 is represented as, 
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with denotes spin and parity of the i th meson and Jj = (0, 
1) is the meson's isospin. Here we neglect the energy trans- 
fer q° carried by the meson for different masses and cut off 
parameters in the real and imaginary parts of the amplitude in 
Eqn. (9). The contribution of i th -meson to the NN-scattering 
amplitude by taking all kinematic is given as, 
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Here the direct and exchange momentum transfer are q = k Q — ko and Q = k l — k\. The first term in Eq.(13), which is already 
of the form of Eq.(9), can easily identify the contribution of T L . The second term is unlike to this form, so we rewrite this as, 
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where the transformation matrix is given as, 
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The row and columns are labeled in the order of S, V, T, A, forms are written as, 
PS. The contribution to the Lorentz invariants (J-l) m simpler 
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Here E c is the total energy in the NN center of mass sys- 
tem. Note that /* depends only on the magnitude of the 
three momentum transfer and the expressions are used to 
fit NN-scattering amplitude at laboratory energy. The full 
parametrizations are frame out in Refs. JUtHtl. 



IV. NUCLEON-NUCLEUS OPTICAL POTENTIAL 

The Dirac optical potential U op t (q, E) can be written as, 

A 



U opt {q,E) = -^(i>\J2^P(k-x(n))T(q,E;n)m20) 



is the nuclear ground state wave function for A-particle, q is 
the momentum transfer and E is the collision energy for a sta- 
tionary target (nucleus) and incident projectile (proton). In the 
present calculation the nuclear recoil energy is neglected be- 
cause of elastic scattering. The operator J-(q, E; n) describe 
the scattering of the projectile from target nucleon 'n' without 
separation into direct and exchange terms. Let us define the 
nuclear ground state by a Dirac-Hartree wave function ll28ll 
and the incident projectile wave function as U{x), then the 
optical potential on incident wave projected to the co-ordinate 
space can be written as, 
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where T is the scattering operator, p is the momentum of the 
projectiles in the nucleon-nucleus center of mass frame, \tp) 
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x\U opt \U ) = / d 3 yd 3 yd 3 xU a (y ) x {(xy \t(E)\x y) + (-if (y x\t(E)\x y)} U (x )U a {y). (21) 



The antisymmetrised matrix element of t(E) in coordinate space is the Fourier transforms 12811 of the matrix element in the 
momentum space co-ordinate and is written as, 
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and similarly for the exchange part tj^(\Q\; E). The nuclear 
density is defined by a simple expression similar to the equa- 
tion of RMF and E-RMF density, 



The first term in the Eqn. (22) defines the direct optical poten- 
tial, 



p L (x,x) = ^U a \ L U a ,p L {x) = p L (x, 



(25) 



The prime stands for occupied states, i.e., sum over target pro- 
tons (pp-amplitude) and target neutrons (pn-amplitude) used. 



U%(r,E) 



-Amp 
M 



d 6 x p L (x )t^{\x~x \;E), (26) 



The nonlocal second term is treated in nonlocal density ap- 
proximation 112911 . which contains plane wave status for inci- 
dent and bound nucleons. We replaced the exchange integral 
with local potential by, 
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where jo is the spherical Bessel-function. The off diagonal one body density is approximated by the local density which result 

as, 
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with s = \x — x | and kf is related to the nuclear baryon 
density by p B (l/2(x + x')) — 2kj/3ir 2 . Now the optical 
potential have the form, 



U opt = U S + j°U v - 2ia.rU 



where 



U L = U L (r,E) = Uk(r,E) + U^(r,E). 



(29) 



(30) 



As the tensor contributions are small, by neglecting these, the 
Dirac equation for projectile has precisely the similar form as 
in RMF and E-RMF equation. By taking the Fourier Trans- 
form of this equation, we get the optical potential as, 
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This equation includes all meson exchanges (except the pseudoscalar meson) with derivative coupling, which is written in the 
form, 



d^q q 2 
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The optical potential is modified by Pauli blocking factor lf30l — 
l34ll a(E) with local density approximation as follows, 
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Here ps is the local baryon density of the target and po is the 

nuclear matter density at saturation. The approximation de- 

2/3 

pends on p B , which agree with phase-space arguments based 
on isotropic scattering. The detail about the Pauli blocking 
factor is given in Ref. I! 811 . To solve the scattering state Dirac 
equation, the wave function is separated into two components 
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(upper and lower) and this equation is expressed as two cou- 
pled first order differential equations. Elimination of the lower 
component leads to a single second order differential equation 
with spin-orbit as well as both local and nonlocal potential. 
The nonlocal Darwin potential can be separated by rewriting 
the upper component of the wave function, A 1/2 (r, E)U(x) 
and 



A(r,E) = 1 



U s (r, E) - U v (r,E) 
E + M ' 



(34) 



After some algebra, the equation can be written as, 
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where the energy-dependent optical potentials are 
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Since the two component Dirac wave functions are eigenstate 
of a.L, so by taking the second derivative of the function we 
can solve easily using Numerov algorithm ll35ll36ll . Note that 
(36) U{x) is not equal to the upper component wave function in the 
region of the potential but when A(r, E) — >l, as r — > oo 

= — i (37) and U has the same asymptotic behavior the wave function 

at large r. Thus the correct boundary condition is imposed 
by matching U to the form of Coulomb scattering solution 
(38) incident in the z-direction 13711 . 
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with E = \/p 2 + M 2 , Xinc is a two-component Pauli spinor, 
6 is the scattering angle, n is the normal to the scattering plane 
and n = Ze 2 /p 2 with Z is the nuclear charge. The scattering 
observables like differential scattering cross-section (^) and 
other quantities, like optical potential (JJ op t), analyzing power 
(A y ) and spin observables (Q— value) are easily determined 
from the scattering amplitude, which are written as, 
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(40) 
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(i) we generate the complex NN-interaction from the Lorentz 
invariant matrix J- (q, E) as defined in Eq. (2). Then the in- 
teraction is folded with the ground state target nuclear density 
for both the RLF (nl [H and MRW parameters OH] sepa- 
rately and obtained the nucleon-nucleus complex optical po- 
tential U op t(q, E). It is to be noted that the pairing interaction 
has been taken into account using the Pauli blocking approx- 
imation. Here, the Pauli blocking enters through the interme- 
diate states of the f-matrix formalism, which has geometrical 
effects on the optical potential, (ii) we solve the wave function 
of the scattering state utilising the optical potential prepared in 
the first step by the well known Numerov algorithm 13511 . The 
result is approximated with the non-relativistic Coulomb scat- 
tering for a wide range of radial component which yields the 
scattering amplitude and other observables lf37ll . By compar- 
ing our calculations with the available experimental data, we 
examine the validity of our RIA predictions for describing 
A y and Q-values which are presented in Figures 1 — 11. 



V. DETAILS OF THE CALCULATIONS AND RESULTS 



First we calculate both the scalar and vector parts of the 
neutrons and protons density distribution for 4 °. 42 > 44 . 48 c a 
from the RMF (NL3) and E-RMF (G2) formalisms Hi]. Then 
evaluate the scattering observables using these densities in the 
RIA frame-work[38], which involves the following two steps: 



A. The neutron and proton densities 

In Fig. 1, we have plotted the proton p p and neutron p n 
density distribution for ^^^^Ca using NL3 and G2 pa- 
rameter sets within RMF and E-RMF formalisms. From the 
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FIG. 1: The proton p p (upper panel) and neutron p n (lower panel) 
density distribution for 40 > 42 . 44 . 48 Ca obtained from RMF (NL3) and 
E-RMF (G2) parameter sets. The experimental [39] p v and deduced 
l40ll p n for 40 ' 48 Ca are also compared. 
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FIG. 3: The relative isotopic neutron density differences Ap n (r) 
for p n ( 42 Ca) - p n ( 40 Ca), p n ( 44 Ca) - p n ( 40 Ca), p n { 4S Ca) - 
p„( 40 Ca), p n ( 44 Ca) - Pn ( 42 Ca), p n ( 4S Ca) - p n { 42 Ca), and 
p n ( 4a Ca) - p n ( 44 Ca). The RMF (NL3) and E-RMF (G2) Ap„(r) 
are compared with the density-matrix-expansion (DME) data [43] 
and the uncertainty deduced neutron difference l42Tl . 
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FIG. 2: The relative isotopic proton density differences Ap p (r) for 
p P ( 42 Ca) - p p ( 40 Ca), p p { 44 Ca) - p p { 40 Ca), and p p ( 4S Ca) - 
p p I 40 Co) obtained from RMF (NL3) and E-RMF (G2) are compared 
with the data [41] in (a), (b) and (c), respectively. 
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FIG. 4: The real (Re) and imaginary (Im) part of the scalar S and 
vector V Dirac optical potential is plotted as a function of the nuclear 
radius for p + 40 Ca system using RMF (NL3) and E-RMF (G2) 
densities, (a) is for RLF and (b) is for MRW parametrization. The 
energy of the projectile proton is Ei a b = 1000 MeV. 



figure, we note that, there is a very small difference in the 
densities for NL3 and G2 parameter sets. However, a care- 
ful inspection shows a small enhancement in central density 
(0-1.6 fm) for NL3 set. On the other hand the densities ob- 
tained from G2 is elongated to a larger distance towards the 
tail region and this nominal difference has significant role to 
play in the scattering phenomena, which is explained later on. 



Further, the agreement of £>„ with the experiment ll39ll and p n 
with the deduced data B4011 for NL3 set is slightly better than 
that of G2. Explicitly, it is worth mentioning that the p p (NL3) 
matches with the data even at the central region, whereas the 
p p of G2 under-estimates through out the density plot. 

A microscopic investigation of Fig. 1 shows a change in 
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FIG. 5: The elastic differential scattering cross-section as a 
function of scattering angle c .m. (deg) for 40 C 'a using both RLF and 
MRW parametrizations at Ei a t = 300 MeV. Triangles are the exper- 
imental data 1 44]. 
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FIG. 6: Same as Fig. 5 for to,*2M,4* Ca at Em _ 80 o MeV. The 
experimental data are taken from |45ll . 



Pp(t), Pn(r), i.e. the area covered by the proton and neu- 
tron densities gradually increases with the mass number in an 
isotopic chain. From the p v (r) and p n (r), we estimate the 
possible relative isotopic density difference Ap(r) for RMF 
(NL3) and E-RMF (G2) parameter sets (see Figs. 2 and 3). 
The calculated Ap p (r) are compared with the experimental 
data Ir4lll in Fig. 2. The measured data of Ap p (r) lies in be- 
tween the prediction of NL3 and G2 values as shown in Fig. 
2. Comparing p P ( 42 Ca) ~ p p { 40 Ca), p P ( 44 Ca) - Pp ( 40 Ca), 
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FIG. 7: Same as Fig. 6 at i?i a f, = 1000 MeV. The experimental data 
are taken from I46fl . 

and p p ( 48 Ca) - p p ( 4Q Ca) of Fig.2 [(a), (b) and (c)], we no- 
tice a better agreement of NL3 values over G2 with respect to 
experimental measurement in the isotopic chain which subse- 
quently reflects in the results of scattering observables. 

The relative isotopic density difference for neutron Ap n (r) 
is compared in Fig. 3 with the deduced neutron density dif- 
ference data 14211 and the density-matrix-expansion prediction 
Il43ll . The predicted results with RMF (NL3) agree well only 
for the double closed shell nuclei 40 C a and 48 Ca. But in 
case of E-RMF (G2) we get excellent match with the deduced 
Ap n (r) for the considered isotopic chain. There is a peak ap- 
pears in Ap n (r) at radial range r ~ 3.4— 3.8 fm and this peak 
slightly shifted towards the center with the increase of neutron 
number. Although Ap n (r) for G2 set gives better agreement 
with the deduced values, the use of NL3 set in the RIA formal- 
ism works well for the scattering observables (shown later). 



B. Optical potential 

With the densities in hand, we calculate the optical potential 
Upot for 4 °> 42 > 44 > 48 £7 a by folding the density matrix with the 
NN scattering amplitude of the proton projectile for 300, 800 
and 1000 MeV. The U po t is a complex function which con- 
stitute both real and imaginary part for both the scalar and the 
vector potentials. In Fig. 4, we present the U po t for ~J? + 40 Ca 
at laboratory energy = 1000 MeV as a representative 
case. We also examine the U po t for other Ca isotopes and 
find similar trends with p + 40 Ca. In other words, we do 
not get any significant difference in the optical potential with 
the increase of neutron number. Similar to the density dis- 
tribution in NL3 and G2 (Fig. 1), here we find a difference 
in U op t(q, E) between the RLF and MRW parametrizations. 
The evaluation methods of the optical potentials using RLF or 
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FIG. 8: Analyzing power A v and Q- Values (spin observables) as a 
function of scattering angle # c .m.(deg) for 40 Ca at Ei a b ~ 300 MeV. 
The experimental data are taken from 14411 . 



FIG. 10: Analyzing power A y as a function of scattering angle 

0c.m.(deg) for 4 0. 4 2,44,48 (7a at = 10 qq Me y 



E, = 800 MeV 

lab 




FIG. 9: Same as FigJS, but at E iab = 800 MeV. The experimental 
data are taken from 14511 . 
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FIG. 11: The spin observable Q— value as a function of scattering 
angle 6> c . m .(deg) for «,42,44,48 Ca at ^ = 1000 MeV 



MRW (see Fig. 4) are somewhat different from each other, 
which is given in Appendix A 13811 . which is responsible for 
the use of the different parametrizations at various ranges of 
incident energies. For example, the RLF parameters used here 
are from Refs. II171 1 1 811 which are computed for energies up 
to 400 MeV and are therefore suitable for lower Ei a i, whereas 
the MRW is better for the higher values which will be dis- 
cussed in the coming sections. Further, the U opt (q, E) values 
from either RLF or MRW, differs significantly depending on 
the NL3 or G2 force parameters. That means, the optical po- 



tential is not only sensitive to RLF or MRW but also to the 
use of NL3 or G2 densities. Investigating the figure, it is clear 
that the extreme values of the magnitude of real and imaginary 
part of the scalar potential are -382.9 and 1 10.6 MeV for RLF 
(NL3) and -372.4 and 177.8 MeV for RLF (G2) respectively. 
The same values for the MRW parametrization are -217.7 and 
40.2 MeV with the NL3 and -333.8 and 61.7 MeV with the 
G2 sets. In case of the vector potential, the extreme values 
for the real and imaginary parts are 293.0 and -136.0 MeV 
for RLF (NL3) and 319.7 and -157.5 MeV for RLF (G2) but 
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with MRW parametrization these appear at 124.1 and -82.3 
MeV with the NL3 and 115.5 and -77.1 MeV with the G2. 
From these variations in the magnitude of scalar and vector 
potentials, it is clear that the predicted results not only depend 
on the input target density, but also sensitive to the kinemat- 
ics of the reaction dynamics. A further analysis of the re- 
sults for the optical potential with RLF, it is noticed that the 
U op t value extends for a larger distance than MRW. For ex- 
ample, with RLF the central part of U op t is more expanded 
than MRW and ended at r ~ 5/m, whereas the U opt persists 
till r ~ 6/to. It is important to point out that the lack of 
the availability of experimental data for optical potential, we 
are unable to justify the capability of parametrizations at dif- 
ferent energies. We also repeat the calculations without Pauli 
blocking and found almost identical results for optical poten- 
tial at Ei ab ~ 300, 800 and 1000 MeV. The effects of RLF and 
MRW parametrizations are presented in the next subsections 
during the discussion of scattering observables. 

C. Differential scattering cross-section 

Evaluation of the differential elastic scattering cross-section 
defined in Eqn. (40) is crucial to study the scattering phe- 
nomena. The results of our calculation for ~f + 40 Ca and 

^40,42,44,48 (j a svs t ems a t incident energies 300, 800 and 
1000 MeV, respectively are displayed in Figs. 5, 6 and 7 along 
with the available experimental data 1144144611 . As it is stated 
earlier that the RIA prediction with the NL3 density is better to 
the choice of G2 for all the angular distributions, irrespective 
of the use of RLF or MRW parametrizations. Again consider- 
ing the energy of the projectile, the RLF predictions best fit to 
the data for Ei a b < 400 MeV (see Fig. 5). However, results 
obtained from the MRW parametrization is better for higher 
incident energies (Figs. 6 and 7) (E lab > 400 MeV) Ojllll]. 
This result shows a fundamental difference between the RLF 
and MRW parametrization depending upon the incident en- 
ergy ranges. Perhaps due to this reason, the explicit off shell 
behavior of RLF and MRW is drastically affecting the scatter- 
ing predictions. Similarly for the optical potential the results 
are insensitive to the Pauli blocking. 



The calculated A y and Q— values obtained by these two forces 
differ significantly from each other for the choice of RLF and 
MRW parametrizations. Also, we observe small oscillations 
in the values of A y and Q with the increase in scattering angle 
9cm. for both RLF and MRW. This oscillatory behavior could 
be related with the dispersion phenomenon of the optical po- 
tential. Similar to the here also the prediction of MRW 
is best fitted to the data for the higher and RLF for lower in- 
cident energies. Further, investigation into the spin rotation 
parameter Q— value, the peak shift and diminished magnitude 
with the increase in neutron number (see Fig. 11) agrees with 
the calculation of first order Brueckner theory using Urbana 
V14 soft core inter-nucleon interactions [47]. It leads to the 
nucleon finite size correction more realistic and hence merits 
a structure effect for the formation of exotic nuclei in labora- 
tory. 

VI. SUMMARY AND CONCLUSION 

We have calculated the density distribution of protons and 
neutrons for 40,42,44,48^ by using jyyrp (NL3 ) and E-RMF 

(G2) parameter sets. From these densities, we estimate the 
relative isotopic neutron density difference for both the force 
parameters. The comparison of Ap n (r) with the data lT42tl 
indicates the superiority of G2 over NL3. The small differ- 
ence in the density at the central region significantly affect 
the results of scattering observables including the optical po- 
tential. A fundamental difference between RLF and MRW 
parametrizations as well as RMF (NL3) and E-RMF (G2) sets 
in the RIA predictions is noticed from the observation of 
A y and Q-value. We conclude from our calculations that RLF 
relatively works well at lower and MRW at higher incident 
energies. The predicting capability of scattering observables 
of RMF (NL3) over E-RMF (G2) is also realised. 

In conclusion, the reaction dynamics highly depends on the 
input density and the choice of parametrization. In addition 
to this, our present study indicates that the RIA is a power- 
ful predictive model which provides a clear picture about the 
successful Dirac optical potentials and can be useful to study 
both stable and exotic nuclei. 



D. Analyzing power and Spin Observable 

The analyzing power A y and the spin observable (Q- Value) 
are calculated from the general formulae given in Eqns. (41) 
and (42) respectively. The results of our calculations for 
p + 40 Ca system at incident energies 300 MeV and 800 MeV 
are shown in Figs. 8 and 9. The RIA predictions for A y using 
RLF with RMF (NL3) density show a quantitative agreement 
with the data [44] at 300 MeV whereas this observation is just 
reverse at 800 MeV lT45ll . That means, the prediction of A y 
resemble the 4£ observations of Figs. 5 — 7. In Figs. 10 and 
11, we present the A y and Q-value for ~f +40,42,44,48 Ca 
composite system at 1000 MeV. These results are obtained 
for both the RLF and MRW parametrizations with NL3 and 
G2 densities in comparision with the experimental data [46]. 
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Appendix A 

If RLF is our choice, the t L functions in Eqns. (17-19) 
and (23-24) involves all the occupied states for pp and pn 
scattering. It is most convenient to shift variables from x' — > 
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x + x' so the t L functions are not involved in the angular can be written as l38ll . 
integration. Now, the first order optical potential Eqn. (20) 



U L (r, E) 



-Atrip 



M 



d 3 r p L (x + x')t%{r , ;E) 



d 3 r'p L (x + x',x)t^(r';E)j (pr') 



(Al) 



after tfi integration, this become 



U L {r, E) 



-8ir 2 ip 



dr'tZ(r';E) J p L (x + x')duj 
dr't%(r';E)j Q (pr') J p L (x + x', x)du> 



(A2) 



where u> — cos9, (\x + x'\ 2 ) = (r 2 + r' 2 + 2uirr') and 
^Jx+xj ^ _ i _|_ ^j rj j _|_ The integral evaluated by 

Gauss-Laguerre quadrature. At the point (x + x'), the radial 
integration must go roughly twice the nuclear radius. Note 
that for spherical nuclei only the scalar and vector are taken 



into account, as the tensor terms are negligible. 

In case of MRW, the optical potential U op t is calculated 
somewhat differently from the RLE Here we tranform the 
density p L (x) to momentum space, then multiply with the 
J rL (q, E), and back which leads to the equation 



U L (r,E) 



M 



(A3) 



with T L {q,E) = F^{E)e- q ? ^ at each proton energy E. The final equation is obtained by adding the contributions from 
proton and neutron states to the direct term Eqn. (A3) which is given as, 



U L (r,E) 



Sip 1 
r 



M 



dqSin(qr) / dr'r' Sin(qr')rf {E)e~ q 13 {E) p L (r') 

u Jo 



(A4) 



r 



The above integrals is solved by double Gussian summation 
methods. 
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